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Abstract 


Signals  from  radio  star  sources  and  satellites  are  regularly  used  for  study¬ 
ing  ionospheric  irregularities.  Amplitude  and  phase  deviations  can  be  imposed 
on  the  signals  from  these  sources  as  they  traverse  the  ionosphere.  A  parameter 
frequently  used  to  describe  the  magnitude  of  this  scintillation  effect  is  scintilla¬ 
tion  index. 

An  experiment  was  designed  to  correlate  various  methods  of  making  scintil¬ 
lation  measurements;  observations  were  made  at  40  MHz  using  the  ionospheric 
beacon,  S-66.  It  is  shown  that  when  the  law  of  the  receiver  detector  is  known, 
a  conversion  method  allows  comparison  of  data  and  statistics  on  scintillation 
index. 

A  simplified  method  of  scaling  scintillation  index  is  described.  The  accu¬ 
racy  of  the  simplified  method  is  determined  by  a  comparison  with  measurements 
of  scintillation  index  by  machine  computation. 
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A  Proposed  Index  for  Measuring 
Ionospheric  Scintillations 


I.  INTRODUCTION 

Studying  the  effects  of  ionospheric  irregularities  often  involves  measuring 
the  amplitude  and  phase  variations  imposed' on  a  signal  as  it  traverses  the  iono¬ 
sphere.  Many  workers  have  used  signals  from  radio  stars  and  satellites  to 
record  the  changes  in. amplitude  known  as  scintillations.  Several  measures  are 
used  to  characterize  the  depth  oPscintillations;  each  of  these  measures  is  called 
a  "scintillation  index"  (S.  I.  ). 

One  measure  is  to  use  a  scale  from  0  to  5  and  by  visual  inspection,  without 
actual  measurement,  assign  a  value  to  the  record.  Other  ways  involve  scaling 
the>deviation  of  the  signal  amplitude  from  the  mean  amplitude.  Four  measures 
are  thus  possible,  depending  on  whether  mean  deviation  or  root -mean -square 
deviation  is  used  and  whether  the  record  is  proportional  to  voltage  or  power.  If 
the  probability  distribution  of  the  amplitude. deviation  is  known,  it  is  possible  to 
relate-The  four  measures  of  scintillation  index  (Briggs  and  Parkin,  1 963).  While 
it  is  difficult  to  relate  these  measures  theoretically  because^the  probability  dis¬ 
tribution  of  the  amplitude  deviation  is  not  in  general  known,  experimental  compar¬ 
isons  show  that  they  are  proportional  to  each  other.  These  relationships  and 
additional  considerations  are  given  in  the  cited  reference. 
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However,  in  most  cases,  only  a  relative  measure  of  scintillation  index  is 
required.  For  a  statistical' analysis  that  involves  a  long  period  of  time,  such  as  a 
study  of  the  seasonal  dependence  of  scintillations,  it  is  important -that  a  standard 
method  of  data  scaling  be  used.  The  notation  that  has  been  adopted  by  AFCRL  and 
the  JSSG  (Joint  Satellite  Studies  Group)  is: 

Scintillation  Index  = 

P  +  P  . 
max  nun 

where  Pmax'i®  the  power  amplitude  of  the  third^peak  down  from  the  maximum 
excursion  of  the  scintillations,  and  Pm-n  is  the  power  amplitude  of  the  third  level 
up  from  the  minimum  excursion.  The  use  of  this  expression  is  demonstrated  in 
Figure  1.  The  chart'  record  is  calibrated  in  dBm  which  is  converted  to  relative 
power  levels  for  calculation  of  S.  I.  Over  the  period  1730  to  1745  E.S.  T.  an 
index  of  40  percentds  obtained. 


1730 


min 


± 


.32—  115 
J , _ 120 


1745  E.S  .T. 


Figure  1.  136  MHz  Scintillations  from  the  Early  Bird  Synchronous  Satellite 

c  ,  .  -max  :  pmia.  _ 

b,1‘  '  Pmax+Pmin"  40 


However,  even  though  only  relative  values  are  usually  required  in  a  study  of 
ionospheric  scintillations,  problems  do  arise  in  the  comparison  of  data  taken  by 
various  workers.  Some  measures  of  S.  I.  may  have  been  scaled  from  chart  divi¬ 
sions  rather  than  by  using  an  accurate  calibration.  The  S,  I,  scaled  from  the  chart 
divisions  under  an  assumption  thatche  chart  display  is  either  linear  in  power  or 
voltage  can  lead  to  errors;  in  reality  the  chart  amplitude  could  be  somewhere  in 


between,  because  it  depends  on  the  law  of -the  receiver  detector.  Also,  statis¬ 
tics  on  satellite  scintillations  may  be  in  the  form  of  a  voltage-index  distribution 
based  on  a  voltage  calibration. 

For  comparison  of  data,  it  is  useful  to  have  a  conversion  relationship  between 
power  index,  S.  I.  and  voltage  index,  S.  I.  If  dat~  was  obtained  from  a 
receiver -with  AGC  and  the  slope  of  the  detector  curve  is  known,  scintillation 
index  values  based  on  chart  divisions  can  be  converted  for  comparison  on  a 
S.  I.  basis. 

Because  of  the  emphasis  on  relating  statistics  on  ionospheric  scintillations 
to  the  problems  of  communication  engineers,  the  following  method  was  developed 
for  converting  data  to  a  common  base,  S.  I.  ^  pj.  The  amount  that  a  signal  fades  is 
important  to  communication  engineers  and  is.  related  to  scintillation  index  by  the 
following; 


Signal  fades  (dB)  =  10  log  1-S.  I. 


=  20  log  jl-S.sI. 


For  example,  if  S.  I.  ^pj  =  0.  5  =  50  percent,  then  the  signal  fade  is  3  dB. 

Since  the  signal  fade  in  dB  would  be  the  same  on  either  a  linear  power  or 
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linear  voltage  display,  the  relationship  is  formed:  S.  I.  ^  =  2  S.  I.  ^  -  S.  I. 

This  conversion  method  can  be  extended  to  cover  receivers  that  are  not 
linear  with  input  power  or  voltage,  such  as  a  receiver  with  AGC,  but  do  have'.a 
constant  slope  to  the  detector  curve  over  the  fading  range,  by  the  following: 


1  —  ( 1  -S.  I.  (v)  ]X 


where  S.  I.  ^  is  the  scintillation  index  measured  in  chart  divisions  with  a  receiver 
that  has  a  detector  law  exponent  of  x. 


2.  KXI’KRIMKNTAI.  RKSUI.TS 

To  determine  the  accuracy  of  the  described  method  for  converting  S.  I.  values 
to  a  con- mon  base,  the  following  experiment  was  made.  Three  R-390  Collins 
receivers  were  gain  adjusted  so  that  a  40  MHz  satellite  signal  received  on  a  com¬ 
mon  antenna  and  a  common  frequency  converter,  whose  output  had  been  equally 
divided  into  three  parts,  would  detect  the  identical  signal  in  three  different  modes: 
linear  power  mode,  linear  voltage  mode,  and  the  AGC  mode.  As. shown  in 
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Figure  2,  the  detector  characteristics  of  each  receiver  had  been  carefully  mea¬ 
sured  to  determine  over  what  signal  range  the  receivers  would  be  essentially 
linear  with  power  or  voltage  or  maintain  a  constant  detector  law. 


By  careful  adjustment  of  the  receiver  gains  and  the  chart  recorder  gains  to 
maintain  a  full-scale  deflection  for  the  same  input  level,  it  was, possible  to 
simultaneously  record  identical  signal  excursions  in  the  three  modes.  Each 
channel  was  accurately  calibrated  so  that  scintillation  index  could>be  measured 
for  the  three  modes  for  identical  signal  excursions  either  using  chart-division 
scaling  or  the  voltage  and  power  calibration.  A  sample  chart  record  is  shown 
in  Figure  3.  For  this  analysis  identical  signal  peaks  were  used  for  the  calcula¬ 
tions  of  S.  I. 
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iFigure  3.  Illustration  of  Scaling  of  Scintillation  Index  for  Three  Detector  Char¬ 
acteristics 


Scintillation  index  was  calculated  for  various  degrees  of  scintillation  for  all 
three  modes.  These  index  numbers  were  then  plotted  against  each  other  in  vari¬ 
ous  ways  to  determine  the  practicability  and  relative  statistical  accuracy  of  com¬ 
paring  (1)  results  of  equipment  having  dissimilar  characteristics,  (2)  scintillation 
measurements  derived  when  calibrations  were  not  taken  and  a  detector  charac¬ 
teristic  law  was.  assumed,  and  (3)  conversion  of  scintillation  measurementsto 
a  common  standard. 

The  results  of  the  experiment  to  compare  scintillation  index  measurement 
are  shown  i.i  Figures  4  through  8.  Figure  4  shows  the  plot  of  the  index  obtained 
by  using  the  signal  power  calibration  for  the  linear  voltage  receiver  (middle  chan¬ 
nel  of  Figure  3)  vs.  the  index  for  the  linear  power  receiver  (lower  channel).  If 
perfect  accuracy  was  obtained  in  the  calibration  and  scaling,  each  channel  would 
indicate  the^same  S.  I, ;  ihowever,- there  is  a  mean  error  for  the  plotted  points  of 
about  +  2.  2  percent. 
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Figure  4.  Measured  Values  of  Scintillation  Index  from 
Linear  Voltage  Receiver  vs.  Measured  Values  of  Scintil¬ 
lation  Index  from  Linear  Power  Receiver.  Power  cali¬ 
bration  for  both  receivers 

Figure  5  shows  the  plot  of  the  index  obtained  by  using  the  signal  power  cali¬ 
bration  for  the  AGC 'receiver  (upper  channel  of  Figure  3)  vs.  the  index  for  the 
linear  power  receiver.  There  is  a  mean  error  for  the  plotted  points  of  +  4  per¬ 
cent. 

Figure  6  shows  the  plot  of  the  index  obtained  by  using  the  signal  voltage 
calibration  for  the  linear  voltage  receiver  vs.  the  index  for  the  linear  power 
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receiver.  The  solid  curve  represents  the  relationship  S.  L  (p,  =  2  S.  I.  (y)  -S.I.(v), 
The  scatter  in  the  points  shows  a  mean  dispersion  of  about  4  percent. 

Figure  7  represents  the  same  data  as  Figure  6,  but  in  this  case  the  voltage 
index  was  scaled  using  chart  divisions. 

Figure  8  shows  the  plot  of  index  values  taken  from  the  AGC  channel  using 
chart  divisions.  The  solid  curve  represents  the  relationship  S.  I.  ^  =  1  - 
( 1  -S.  I.  (agc:))4'  5*  T**e  exPonent  5  represents  the.mean  slope  of  the  AGC  - 
detector  characteristic  as  shown  in  Figure  2. 
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Figure  5.  Measured  Values  of  Scintillation  Index  from  Receiver  with  AGC  Char¬ 
acteristic  vs.  Measured  Values  of  Scintillation  Index  from  Linear  Power  Receiver. 
Power  calibration  for  Doth  receivers 
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Figure  6.  Measured  Values  of  Scintillation  Index  from  Linear  Voltage  Receiver 
vs.  Measured  Values  of  Scintillation  Index  from  Linear  Power  Receiver.  Volt¬ 
age  calibration  for  linear  voltage  receiver.  Power  calibration  for  linear  power 
receiver 
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Figure  7.  Measured  Values  of  Scintillation  Index  from  Linear  Voltage  Receiver 
vs.  Measured  Values  of  Scintillation  Index  from  Linear.  Power  Receiver.  Chart 
division  calibration  for  linear  voltage  receiver.  Power  calibration  for  linear 
power  receiver 


Figure  8.  Measured  Values  of  Scintillation  Index  from  Receiver  with  AGC 
Characteristic  vs.  Measured  Values  of  Scintillation  Index  from  Linear  Power 
Receiver.  Chart  division  calibration  for  AGC  receiver.  Power  calibration 
for  linear  power  receiver 
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3.  scai.im;  n-*  scintii iniiex 

It  was  pointed  out  in  a  preceding  section  that  scintillation  index  was  defined  as 

P1112*  +  p”?-n  -  At  first  the  measure  of  S.  I.  was  made  by  measuring  the  two 
max  min 

levels  and  then  computing,  either  fay  machine  or  fay  hand,  the  resulting  S.  I. 

The  first  step  was  to  mark  the  P  and  P  .  levels  on  the  chart  for 

max  min 

each  time  period  that  S.  I.  was  to  be  read.  P  was  arbitrarily  chosen  as  the 

max 

third  peak  dow'n  from  the  maximum  and  P  .  as  the  third  minimum  up  from  the 

min 

deepest  fade.  In  the  case  of  the  low  altitude  satellite,  S.I.  was  read  either  at 
the  peak  of  the  Faraday  cycle  or  once  per  minute  if  Faraday  periods  were  ob¬ 
scured  by  heavy  scintillation.  For  the  synchronous  satellite  signals,  scintilla¬ 
tion  index  was  measured  in  each  5  minute  interval.  The  time  constant  of  the 
recording  system  -was  chosen  to  be  short  compared  with  the  fastest  scintillation. 

The  second  step  was  to  read  P  and  P  .  as  a  relative  power  level  based 

max  mm 

on  an  accurate  amplitude  calibration  of  the  chart  deflection.  If  the  calibration  had 

been  made  in  decibel  steps,  it  was  converted  to  an  appropriate  numerical  ratio. 

Scintillation  index  was  either  computed  manually,  using  the  formula,  or  Pmax 

and  P  .  were  punched  on  cards  for  machine  calculation  of  S.I.  and  further 
mxn 

analysis. 

•A  simplified  method  of  measuring  S.  I.  has  evolved  which  is  now  implemented 
for  all  S.  I.  measurements.  It  again  depends  on  an  accurate  amplitude  calibra¬ 
tion  in  decibel  steps  of  the  chart  deflection.  P  .and  P  .  are  determined  in 

max  mxn 

the  same  manner  as  described  in  step  1  above.  The  decibel  change  of  P 

b  max 

Pmin  *s  t*len  reac*.  using  the  receiver  amplitude  calibration  curve.  The  mea¬ 
sured  values  in  decibels  of  P  -  P  are  then  converted  to  S.  I. ,  using  the 
curve  shown  in  Figure  9.  If  further  analysis  is  required,  then  either  the  decibel 
values  of  Pmax  -  pmjn  or  the  values  converted  to  scintillation  index  can  be 
punched  on  cards  for  machine  processing. 

The  conversion  graph  shown  in  Figure  9  was  determined  by  assuming  equal 

percentage  changes  from  the  average  level  for  P  and  P  .  ,  changing  those 

max  mm 

values  to  a  decibel  change.and  then  summing  them  for  the  total  decibel  change 

as- shown  in  Table  1.  For  example,  based  on  an  arbitrary  average  level  of  1.  0, 

a  S.  I.  =  50  percent  corresponds  to  a  P  .  =  0.  5  =  3  dB,  aP  =  1.5  = 

min  max 

1.  77  dB,  and*  therefore,  a  total  change  for  P  -  P  .  of  4.  77  dB. 

b  max  mm 

As  a  check  on  the  accuracy  of  the  simplified  method  compared  with  the 
previously  used  method  of  machine  computation,  several  values  of  scintillation 
index  were  measured  using  both  procedures.  The  calibration  of  an  S-66  record 
at  41  MHz  was  chosen  because  the  average  level  for  low  altitude  satellites 
changes  over  a  far  greater  range  than  that  from  synchronous  satellites. 
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The  results  of  this  comparison  are  shown  in  Figures  10,  11,  and  12.  The 

measurements  of  P  -  P  .  in  dB  and  their  conversions  to  S.  I.  are  plotted 
max  mm 

against  the  normal  method  of  machine  computation  of  S.  I.  Figure  10  corresponds 

to  the  case  of  measuring  S.  I.  for  weak  signals,  that  is,  where  the  average  * 

level  is  about  1/4-of  the  chart  width.  Figure  11  for  average  levels  of  about  1/2  of  j 

the  chart  width,  and  Figure  12  for  strong  signals  with  an  average  level  about  3/4  j 

of  the  chart: width.  The  deviations  of  the  points  from  the  dotted  line  show  the 

lack  of  perfect  agreement  between  th4  two  methods  of  measurement.  The  greatest  j 

differences  are  for  the  small  signal  case  and  are  due  to  the  compression  of  the  j 

calibration  curve  at  small  signal  levels  which  results  in  greater  inaccuracies  of 
reading- the  value  of  the  signal. 


Figure  10.  Comparison  of  Scintillation  Index  Measurements.  S.  I.  measured 
directly  from  calibration  and  conversion  curve  vs.  S..I.  measured  with  machine 
computation  of  Pmax  -  Pmin/Pmax  +  Pmin!  small  signal  case— average  level 
approximately  1/4  chart  width  (10-15  mm) 
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S.I.HU  NORMAL  METHOD  -  MACHINE  COMPUTATION 


Figure  11.  Comparison  of  Scintillation  Index  Measurements.  S.  I.  measuied 
directly  from  calibration  and  conversion  curve  vs.  S.  I.  measured  with<machine 
computation  of  Pmax  -  Pmin/Pmax  +  Pmin!  medium  signal  case-average  level 
approximately  1/2  chart  width  (20-30  mm) 


S.I.  (%)  NORMAL  METHOD  -  MACHINE  COMPUTATION 


Fi’gyre  12.  Comparison  of  Scintillation  Index  Measurements.  S.  I.  measured 
directly  from  calibration  and  conversion  curve  vs.  S.  I.  measured  with  machine 
computation  of  Pmax  "  Pmin/Pmax  +  Pminl  large  signal  case— average  level 
approximately  3/4  chart  width  (30-35  mm) 
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If  the  simplified  method  is  to  be  applied  for  analysis  where  there  is  a  small 
signal -to -noise  ratio  and  the  deflection  on  the  chart  due  to  the  sky  background 
temperature  is  appreciable  compared  to  the  signal,  the  receiver  calibration  must 
be  correctly  applied  so  that  only ‘ihe-deflections  from  the  signal  are  being  mea¬ 
sured. 


I.  CONCLUSION 

In  summary,  it  is  felt  that  only  a  relative  measure  of  scintillation  index  is 
necessary  to  describe  amplitude  fluctuations  caused  by  ionospheric  irregularities, 
and  that  the  simple  measure  S.  I.  =  which  was  defined  earlier, 

should  be  used.  To  obtain  the  best  accuracy  when  comparing  scintillation  data, 
the  records  should  be  accurately  calibrated  on  a  power  base.  If  scintillation 
index  is  scaled  from  chart  divisions  and  the-law  of  the  detector  is  known,  the  in¬ 
dex  values  can  be  converted  to  a  power  base  with  a  small  decrease  in  accuracy, 
as  compared  to  scaling  from  a  signal  generator  calibration.  Comparison  of  data 
without  the  benefit  of  a  calibration  or  knowledge  of  the  detector  characteristic 
can  lead  to  errors  of  a  factor  of  2  or  more.  This  is  evident  from  Figure  3 
where  the  apparent  S.  I. ,  as  indicated  by  chart  divisions,  is  very  small  for  the 
AGC  mode  and  la.gest  for  the  linear  power  mode.  The  simplified  method  of 
scaling  S.  I.  gives  adequate  accuracy  for  statistical  studies. 
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